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ENERGY DISTRIBUTION FUNCTIONS OF KILOVOLT IONS PARALLEL 
AND PERPENDICULAR TO THE MAGNETIC FIELD OF 
A MODIFIED PENNING DISCHARGE 
by J. Reece Roth 
Lewis Research Center 

SUMMARY 

The distribution function of ion energy parallel to the magnetic field of a modified 
Penning discharge has been measured with a retarding potential energy analyzer . These 
ions escaped through one of the throats of the magnetic mirror geometry. Simultaneous 
measurements of the ion energy distribution function perpendicular to the magnetic field 
have been made with a charge- exchange neutral detector. The ion energy distribution 
functions are approximately Maxwellian, and the parallel and perpendicular kinetic tem- 
peratures are equal within experimental error. These results suggest that turbulent 
processes previously observed in this discharge Maxwellianize the velocity distribution 
along a radius in velocity space and cause an isotropic energy distribution. The kinetic 
temperatures observed are of the order of kilovolts, and the tails of the ion energy dis- 
tribution functions are Maxwellian for up to a factor of 7 e-folds in energy. When the 
distributions depart from Maxwellian, they are enhanced above the Maxwellian tail. 
Above densities of about 10 10 particles per cubic centimeter, this enhancement appears 
to be the result of a second, higher temperature Maxwellian distribution. At these high 
particle energies, only the ions perpendicular to the magnetic field lines were investi- 
gated. 


INTRODUCTION 

Previous investigations of a modified Penning discharge (ref. 1) revealed ions with 
approximately Maxwellian velocity distributions and kilovolt kinetic temperatures. The 
diagnostic instrument used to measure the energy distribution of ions parallel to the 
magnetic field lines was a retarding potential analyzer, located outside one of the mag- 


netic mirrors. The data reduction procedures and experimental results from an exten- 
sive series of measurements with the retarding potential energy analyzer have been 
reported elsewhere (ref. 2). 

In view of the possible application of the modified Penning discharge to a steady- 
state plasma heating scheme for fusion research, it is desirable to confirm independent- 
ly the high ion kinetic temperatures and Maxwellian velocity distribution observed in the 
ions lost along the magnetic-field lines. This report describes an investigation in which 
the parallel and perpendicular ion energy distribution functions were simultaneously 
measured for 163 experimental runs taken under a range of operating conditions. Such 
simultaneous measurements have not been reported in the literature of Penning dis- 
charges, although Konstantinov, Myskin, Sorokin, and Tsel'nik (ref. 3) have measured 
the energy of charge- exchange neutrals from the midplane of a modified Penning dis- 
charge. In the investigation described in this report, a charge -exchange neutral detec- 
tor was used to measure the perpendicular ion energy, while simultaneous observations 
of the parallel ion energy were made with the retarding potential energy analyzer. 


Characteristics of Modified Penning Discharge 

The operating characteristics of the discharge have been described in references 
1 and 4 ; a detailed description of the superconducting magnet facility in which the dis- 
charge is operated is described in reference 5; and some of the physical processes re- 
sponsible for ion heating are discussed in reference 6. An isometric cutaway drawing of 
the experimental configuration is shown in figure 1. The plasma is approximately 
15 centimeters in diameter at the midplane. The magnetic field consists of a magnetic 
mirror with a 2. 62:1 mirror ratio generated by two superconducting coils, the Dewars 
of which are approximately 18 centimeters in diameter at the magnetic mirror throats. 
The magnet system shown in figure 1 is located in a vacuum tank 1 meter in diameter 
and approximately 2 meters long and gives very good experimental and visual access to 
the experimental volume. The parallel Dewars on either side of the experimental vol- 
ume are at liquid- nitrogen temperature. The anode ring is operated at positive poten- 
tials up to 35 kilovolts with respect to the walls of the vacuum tank, and the electrical 
circuit to the power supply is completed by ions which impinge on the grounded tank 
walls. 

The experimental investigation consisted of 163 runs taken with deuterium gas. The 
maximum value of the magnetic field was 1 tesla at the magnetic mirror throats and 
0. 38 tesla at the midplane, in the vicinity of the anode rings. From previous measure- 
ments taken with a Langmuir probe at the throat of the mirror (ref. 6), the electron 
temperature is known to be approximately 50 eV, within a factor of 2, in this plasma. 
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Figure 1. - Isometric cutaway drawing of modified Penning discharge and superconducting magnetic mirror facility. 
Approximate locations of retarding potential energy analyzer and of fast neutral energy analyzer are shown 
schematically. 
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This is as much as an order of magnitude below the ion temperature at the same operat- 
ing conditions. Any atomic deuterium generated in the plasma volume has ample oppor- 
tunity to recombine on the walls before it can return to the experimental volume, as a 
result of the long mean free paths for binary collision processes. Because of the rather 
large vacuum-tank wall area and the cold liquid- nitrogen surfaces on either side of the 
plasma volume, it is unlikely that a significant population of atomic deuterium can exist 
in the plasma. For this reason it was assumed that the neutral background gas consisted 
of molecular deuterium gas. The principal ionized species within the plasma was taken 
to be ionized deuterium atoms D + . It is believed that D + dominated the plasma be- 
cause investigation of the spectrum of electrostatic potential fluctuations yielded a peak 
at the ion cyclotron frequency of D + (when such a peak could be observed), and peaks 
were absent or much fainter by comparison at the ion cyclotron frequency appropriate 
to ionized molecular deuterium Dg. 

The 163 experimental runs spanned a range of the positive anode voltage from 5 to 

-5 -5 

35 kilovolts and neutral background pressures p Q from 3.9x10 to 12.0x10 torr, 

where the neutral deuterium pressure has been corrected for the gage factor of this 

apparatus. It was estimated from previous work (ref. 6) that the ion density at the mid- 
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plane of the modified Penning discharge ranged from approximately 10 to 6x10 par- 
ticles per cubic centimeter in this investigation. 


Diagnostic Instruments 

The ion energy distribution function was measured by two nonperturbing diagnostic 
techniques, one parallel to the magnetic-field lines, and the other perpendicular to the 
magnetic- field lines. In figure 2(a) is a schematic drawing of the retarding potential 
energy analyzer used to measure the parallel distribution function of ions leaving the 
plasma through the magnetic mirror throat. The analyzer was located approximately 
20 centimeters outside the magnetic mirror throat and 10 centimeters from the magnetic 
axis. The analyzer was parallel to the local magnetic- fie Id line and measured the paral- 
lel component of the ion energy. The raw data, taken on an X, Y-recorder, represented 
an integrated ion energy distribution function. A computer program determined, by 
iteration, the integrated Maxwellian distribution which best fit the raw data, in the sense 
of minimizing the mean square error between the logarithms of the best-fitting distri- 
bution and of the raw data. Using the logarithm of the value of the data point rather than 
the actual value of the data point tended to give approximately equal weight to all portions 
of the data curves. By comparing the best-fitting computer solution to the raw data, 
one may assess the extent to which the raw data are Maxwellian by observing the extent 
to which the best-fitting curve is also a good fit. The principal source of error for this 
application of the retarding potential energy analyzer is in the iterative procedure used 
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X, Y-recorder 


(a) Retarding potential energy analyzer, located approximately 4 centimeters from 
plasma axis and 20 centimeters outside magnetic mirror throat. 
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(b) Schematic drawing of charge-exchange neutral detector. Line of sight of the detector 
was tangent to sheath between plasma and anode ring approximately 6 centimeters below 
magnetic axis. 


Figure 2. - Apparatus for measurement of parallel and perpendicular ion energy 
distribution functions. 





to reduce the data, which sometimes does not converge to the proper kinetic tempera- 
ture. (This was determined by obtaining a best fit to "pseudodata, " an integrated Max- 
wellian of specified kinetic temperature. ) These convergence problems give rise to 
errors in the parallel- ion kinetic temperature which may be as high as 20 percent. The 
analytical theory, error analysis, and experimental application of the retarding poten- 
tial energy analyzer are discussed in reference 2. 

The charge -exchange neutral detector used in these investigations is a geometric 
duplicate of that developed at Fontenay-aux- Roses (ref. 7). The particle detector elec- 
tronic system was developed at the NASA Lewis Research Center. A schematic showing 
this instrument and its location with respect to the plasma is presented in figure 2(b). 
The retarding potential energy analyzer and the neutral detector were positioned to 
sample the anode sheath, approximately 1 centimeter from the inner circumference of 
the anode ring. Fast ions in the sheath between the plasma and the anode ring charge 
exchange on the background of neutral molecular deuterium gas, and the resulting fast 
deuterium atoms then move along straight lines to the walls of the vacuum vessel. In 
the neutral detector, some of the fast neutrals pass through a series of slits and are 
reionized by a nitrogen gas cell, which is maintained at pressures of approximately 
0. 10 torr. The reionized fast particles then pass through additional focusing slits into 
a set of 90° electrostatic deflector plates. The energy of the ions which may pass 
through the 90° deflector plates is determined by the potential applied across them. By 
changing the voltage applied between the deflector plates, one may sweep out the energy 
distribution function of the reionized particles. The ions are passed through a second 
series of slits and impinge on a photomultiplier detector, which counts individual par- 
ticles. 

The particle flux at the photomultiplier is only indirectly related to the ion distribu- 
tion function inside the confined plasma. To obtain the original ion distribution function, 
one must correct for the charge-exchange process of atomic deuterium ions on molecu- 
lar deuterium, 


D + + D 2 - Dg + D° 

One must also correct for the reionization cross section in the nitrogen gas cell and for 
the energy resolution of the electrostatic analyzer. The flux I(V) of individual particles 
reaching the photomultiplier (see the appendix) may be written 

I(V) = C o n i n 0 a x (V)v i f(V)k(V)G a (V) (1) 

where V is the ion energy in electron volts, C Q is a geometrical and dimensional con- 
stant, f(V) is the desired ion energy distribution, k(V) is the reionization probability of 
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atomic deuterium in a nitrogen gas cell, and G a (V) is a correction for the energy reso- 
lution of the analyzer and results from the fact that the relative energy resolution dV/V 
is independent of energy. Thus, the energy resolution involves correction by a factor 

G a (V) ~ V (2) 

The factor n i n 0 ° r x (V)v i is the number of charge-exchange reactions per cubic centime- 
ter per second in the volume sampled by the analyzer, n^ and n Q are the ion and neu- 
tral number density, respectively, and v i is the deuterium ion velocity. The charge- 
exchange cross section a x (V) is available in reference 8 and is known from 70 eV ^ V 
:< 200 keV. The correction factor k(V) for the nitrogen gas cell is known for energies 
over the range 2 < V ^ 100 keV and is 

k(V) ~ V 0, 56 (3) 


over this range (refs. 7 and 8). 

One may convert the raw data for I(V) to the distribution function of ion energy in 
the plasma f(V) by using equations (1) to (3) to obtain 

CjKV) 

f(V) = i (4) 

V 2 - 06 o x (V) 

where is a constant for the apparatus. The ion energy distribution function ob- 
tained from equation (4) was used as input to a computer program which obtained a 
Maxwellian distribution with a least- squares best fit to the experimentally determined 
ion energy distribution f(V). Because the correction factor k(V) was not known for 
nitrogen gas below 2 kilovolts, raw data below 2 kilovolts were not utilized in obtaining 
the best-fitting Maxwellian distribution. The principal error in applying the neutral 
detector arose from establishing the zero level and rounding off errors in reading the 
raw data. These errors, combined with uncertainties in the charge- exchange cross 
sections, may have introduced an error of as much as 10 percent in determining the 
perpendicular ion kinetic temperature. 


EXPERIMENTAL RESULTS 

In figure 3 is shown a characteristic example of data which were taken simultane- 
ously with the retarding potential energy analyzer and the charge- exchange neutral de- 
tector. The anode voltage was 8 kilovolts, and the neutral background gas pressure was 
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(a) Data measured by charge-exchange neutral detector. 
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(b) Data measured by retarding potential energy analyzer. 
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(a) Data measured by charge-exchange neutral detector. 


o Raw data 

A Best-fitting Maxwellian 
V b - 2200 \ln with V|i » 1705 eV 
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(b) Data measured by retarding potential energy analyzer. 


Figure 3. - Example of parallel and perpendicular 
Maxwellian ion energy distribution functions 
measured with instruments shown in figure 2. 
Anode voltage, 8 kilovoltSj deuterium gas pressure, 
5. 2xl0' 5 torr. 


Figure 4. - Example of parallel and perpendicular 
Maxwellian ion energy distribution functions 
measured with instruments shown in figure 2. 
Anode voltage, 35 kilovolts; deuterium gas pressure, 
4. 3x10-5 torr. 


5. 2x10” 5 torr of deuterium. Figure 3(a) is the energy distribution function measured 
with the charge- exchange neutral detector, corrected with equation (4), and plotted as a 
function of energy on a semilogarithmic graph. The corrected raw data are indicated 
by the round symbols, and the Maxwellian distribution which best fits these raw data is 
indicated by the triangular symbols. In figure 3(b) is shown the integrated energy dis- 
tribution function measured with the retarding potential energy analyzer. The round 
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symbols are the corrected raw data, and the triangular symbols are the best- fitting 
integrated Maxwellian distribution. Figure 3(a) gives a perpendicular component of the 
ion kinetic temperature, V ± = 620 eV, and figure 3(b) a best-fitting parallel kinetic tem- 
perature, V (1 = 559 eV. It is characteristic in figure 3 that both the parallel and perpen- 
dicular ion energy distribution functions are approximately Maxwellian, that the Max- 
wellian tail of the neutral detector data extends over approximately seven e- folding 
lengths in energy, and that the parallel and perpendicular ion energies are approximate- 
ly equal within the scatter in the data of the two diagnostic methods employed. 

In figure 4 is a second example of the perpendicular and parallel ion energy spectra, 
under conditions which yielded kinetic temperatures of nearly 2 kilovolts. The operating 
conditions were an anode voltage of 35 kilovolts, and a neutral background gas pressure 

_ R 

of 4. 3x10 torr of deuterium. In figure 4(a) is the neutral detector spectrum with a 
best- fitting perpendicular kinetic temperature of V ± = 1942 eV. In figure 4(b) is the 
retarding potential data, with a best-fitting parallel kinetic temperature of V (( = 1705 eV. 
The knee in the latter data indicates that the magnetic- fie Id line on which the retarding 
potential analyzer was located was floating at a positive potential of V f = 2200 volts, 
an unusually high value for this experiment, but only 6 percent of the applied anode 
voltage. Figures 3 and 4 are characteristic in that it was possible, because of instru- 
ment limitations, to take valid data at higher energies with the neutral detector than 
with the retarding potential energy analyzer. 

Ion energy distributions analogous to those shown in figures 3 and 4 were obtained 
for all 163 runs. Temperatures obtained from the best- fitting Maxwellian distributions 
are shown in figure 5. On the small number of occasions for which a two- temperature 
Maxwellian distribution was observed with the neutral detector, the lower of the two 
kinetic temperatures was plotted in figure 5 . On the abscissa is the kinetic temperature 
derived from the parallel component of ion energy, measured with the retarding poten- 
tial energy analyzer. On the ordinate is the kinetic temperature derived from the per- 
pendicular component of ion energy measured with the charge- exchange neutral detector. 
There is approximate agreement of the parallel and perpendicular components of the ion 
kinetic temperature over the range of data taken. 

The agreement of with and the degree of data spread may be estimated by 
fitting the data in figure 5 to the relation 


V l = C 2 V ll 


(5) 


and obtaining the median and relative standard deviation of C 2 . In figure 6, the cumu- 
lative probability of C 2 for 163 runs is plotted on probability coordinates. The author's 
best estimate of the Gaussian distribution to the data is indicated by the straight line. 



Figure 5. - Paralles and perpendicular ion energies 
determined from charge-exchange neutral detector 
and retarding potential energy analyzer for 163 
experimental runs. 



Proportionality constant, C 2 

Figure 6. - Cumulative probability of observing value of 
proportionality constant C 2 from equation (5) for 
population of 163 experimental runs. Median value of 
C 2 , 1.04; standard deviation, 0.270; relative standard 
deviation, 0.26. 


The median value of C 2 , C 2 , is the intersection of the best-fitting straight line with the 
50 percent line, 

= 1. 04 (6) 

which implies that the perpendicular energy is slightly higher than the parallel energy. 
The relative standard deviation of the data population, the horizontal distance between 
the intersection of the best- fitting straight line with the median and the one standard 
deviation lines, was approximately 


s(c 2 ) 

AC 9 = —~r- - 0.26 (7) 

c 2 

Thus, the parallel and perpendicular ion kinetic temperatures are equal to within the 
estimated experimental error. The parameter C 2 was plotted as a function of the in- 
dependent variables (anode voltage and neutral gas pressure) for the 163 runs, and no 
systematic trend of C 2 was observed. The same physical process which is responsible 
for making C 2 a fixed constant apparently operates over the entire range of parameters 
investigated. 

Two general types of anomalous ion energy distributions were observed in the per- 
pendicular ion energy distribution function. Both were characterized by an increase of 
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the observed distribution function above the best- fitting Maxwellian value in the tail of 
the distribution. This enhancement of the Maxwellian tail assumed two forms. One of 
these was a two- temperature Maxwellian distribution, in which the Maxwellian tail was 
dominated by a hotter component of lower density that gave the type of ion energy dis- 
tribution function illustrated in figure 7. Two- temperature distributions were observed 
only at plasma densities above approximately 10^ particles per cubic centimeter. This 
two- temperature Maxwellian distribution apparently results from the analyzer sampling 
two separate regions with different ion kinetic temperatures. 
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Figure 7. - Example of two-temperature Maxwellian 
distribution observed with neutral detector at ion 
number densities above about 10 3 ® per cubic 
centimeter. Anode voltage, 10 kilovolts; back- 
ground pressure, 7.8xl0'5torr of deuterium gas. 


Two- temperature Maxwellian distributions may also occur as a result of a popula- 
tion of hot impurities superimposed on the distribution function of the principal species 
Such impurities, with each species having a distinct ion temperature, have been de- 
scribed by Stirling (ref. 9). This explanation' is unsatisfactory for the present experi- 
ment, however, since observation of the Penning discharge plasma with a bench spec- 
trometer revealed no detectable impurity lines. Ion energy distributions similar to 
figure 7 have been reported by Artsimovich (ref. 10) in connection with the Tokamak 

series of experiments. Artsimovich reported that two- temperature Maxwellian distri- 
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butions occurred in the Tokamak at number densities below 10 particles per cubic 
centimeter. 
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A second anomalous departure from the Maxwellian distribution was observed in 
approximately a dozen cases, at higher neutral gas pressures. This departure is illus- 
trated in figure 8 and shows an ion energy distribution function with a bump on the Max- 
wellian tail at an energy which roughly corresponds to the voltage applied to the anode 
ring. Such a distribution function may be due to physical processes in the weak plasma 
penumbra which surrounds the modified Penning discharge. 
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Figure 8. - Example of ion energy distribution function, 
measured with neutral detector, with a bump on the 
Maxwellian tail. Anode voltage, 11 kilovolts-, back- 
ground pressure, 1. 2x10'^ torr of deuterium gas. 


DISCUSSION OF RESULTS 

The present investigation has demonstrated generally good agreement between the 
parallel and perpendicular kinetic temperatures observed with nonperturbing diagnostics 
in a modified Penning discharge. These observations provide an independent confirma- 
tion of the results reported in reference 2, in which a sophisticated data reduction pro- 
cedure was applied to raw data obtained from the modified Penning discharge with a re- 
tarding potential energy analyzer. The Maxwellian ion energy distributions observed in 
the distribution functions parallel and perpendicular to the magnetic field are consistent 
with a model in which the turbulent processes described in reference 4 tend to Maxwel- 
lianize the distribution function along a radius in velocity space in such a way that the ion 
energy is also isotropic in velocity space. 

The modified Penning discharge has been shown to be capable of generating a 
steady- state plasma of kilovolt ion energies with an energy distribution function that is 
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Maxwellian and virtually isotropic in velocity space. The modified Penning discharge 
may therefore be of interest in controlled fusion experiments, particularly in those 
toroidal systems in which one desires to create an approximately isotropic distribution 
of hot ions. 

R. P. G. Valckx and A. Bariaud of Fontenay-aux-Roses provided useful conversa- 
tions and assistance with the construction of the fast neutral energy analyzer. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 3, 1972, 

503-10. 
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APPENDIX - SYMBOLS 


B magnetic- field strength, T 

C Q constant appearing in eq. (1) 

Cj constant appearing in eq. (4) 

C 2 constant appearing in eq. (5) 

C 2 median value of C 2 

AC 2 defined by eq. (7) 

6C 2 standard deviation of C 2 

tf T electric field strength, V/m 

f(V) distribution function of ions in plasma 
Ga (V) correction for energy resolution of analyzer 
I(V) particle flux 

k(V) reionization probability in neutral gas cell 

O 

m ion number density, particles/m 

3 

n Q neutral particle density, particles/m 

V ion energy, eV 

perpendicular component of ion energy 
V|| parallel component of ion energy 


v^ ion velocity in plasma 

cr (V) charge- exchange cross section in plasma 

A 
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